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Abstract: We calculate the long-distance effect generated by the four-quark operators 
with c-quarks in the B — t- K^*^i~^i~ decays. At the lepton-pair invariant masses far below 
the cc-threshold, <C 4m^, we use OPE near the light-cone. The nonfactorizable soft- 
gluon emission from c-quarks is cast in the form of a nonlocal effective operator. The B — )• 
K^*^ matrix elements of this operator are calculated from the QCD light-cone sum rules 
with the S-meson distribution amplitudes. As a byproduct, we also predict the charm-loop 
contribution to B ^ K*j beyond the local-operator approximation. To describe the charm- 
loop effect at large q'^, we employ the hadronic dispersion relation with ip = J/ip, i/j{2S), ... 
contributions, where the measured B — )■ K^*^'ip amplitudes are used as inputs. Matching 
this relation to the result of QCD calculation reveals a destructive interference between 
the J/ip and ip{2S) contributions. The resulting charm-loop effect is represented as a g^- 
dependent correction ACg{q'^) to the Wilson coefficient Cg. Within uncertainties of our 
calculation, at below the charmonium region the predicted ratio ACg{q'^)/Cg is < 5% 
for B — )• K£~^£~, but can reach as much as 20% for B — )• K*£'^£~, the difference being 
mainly caused by the soft-gluon contribution. 
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1. Introduction 

Being very prospective channels for the search for new physics, the exclusive B K^*^l'^£^ 
decays have quite a complicated dynamics. In the Standard Model the underlying flavour- 
changing b — )• s£~^£~ transition is described by the effective Hamiltonian [|l|, ^ : 

H,jf = -^vt,v*Y.^,{t,m,,) , (1.1) 

^ 1=1 

a superposition of the effective operators Oj, weighted by their Wilson coefficients Cj and 
normalized at the scale fi ~ rrif,. The relevant operators are presented in Appendix A. In 
this paper we neglect the CKM-suppressed contributions proportional to VubVus^ adopting 
the approximation |14feV;*| ~ \ VcbV*,\ = 0.0428t|]:[}[j^ 
The decay amplitude 

A{B K^*h+£-) = -{K^*h+£- I Heff \ B) , (1.2) 

contains a rich variety of hadronic matrix elements of the operators Oj. The dominant 
contributions to ( |1.2| ) are generated by the 09,10 and Oy^ with large Wilson coefficients. 
In Appendix B the hadronic matrix elements of these operators are presented. They are 
factorized in terms of the B — )• K^*^ form factors, which, similar to the B ^ tt, p form 
factors of the weak semileptonic decays, are obtained from lattice QCD or QCD light- 
cone sum rules (LCSR). In addition, there are specific contributions to the B — )■ K^*^£^£~ 
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and B — )• K*^ amplitudes generated by the current-current operators, Oi^2 and penguin 
operators 03-6,835 combined with the electromagnetic (e.m.) interaction of quarks. A 
major challenge for the theory is to identify and estimate these hadronic matrix elements 
one by one. 

An important effect, which is the main topic of our study, is generated by the current- 
current operators Oi,2 acting together with the c-quark e.m. current. This mechanism 
involves an intermediate "charm-loop" , coupled to the lepton pair via the virtual photon. 
The analogous effects for the u-quark current-current and quark-penguin operators with 
different flavours are not important in i? — )• K^*^£^i~, being either CKM suppressed, or 
multiplied by a small Wilson coefficient. 

The c-quark loop turns into a genuine long-distance hadronic effect if the lepton- 
pair invariant mass = {p£+ + Pe-)"^ reaches the region of charmonium resonances ip = 
{J/ip,ip{2S), ....}. At q"^ = 171"^, the process B — )• K^*H~^£~ transforms into a nonleptonic 
weak decay B — )■ 'il)K^*\ followed by the leptonic annihilation of V'- All i/'-resonances with 

< {mB—'m^(t)) contribute to this mechanism. Moreover, above the DD threshold there 
are continuum contributions of the intermediate charmed hadrons and it is very difficult 
to include them in a model-independent way. Although the (7^-intervals around J/ip and 
V'(2<S') are subtracted from the measured lepton-pair mass distributions in B K^*k+£-, 
the intermediate and/or virtual cc states contribute outside the resonance region and their 
effect has to be accurately estimated. 

Usually, in the leading-order, the c-quark loop diagram is included into the factoriza- 
tion formula for B — )■ K^*^l'^i^ . In addition, hard-gluon exchanges between the c-quark 
loop and the rest of the diagram are taken into account, together with other perturbative 
nonfactorizable effects (see e.g., Q). One generally predicts these effects to be small, if q"^ 
is far below the charmonium region. 

Two important questions concerning the charm-loop effect remain unanswered, despite 
many dedicated studies. The first question is: how important are the soft gluons emitted 
from the c-quark loop and violating the factorization? The second question concerns the 
validity of the approximation "c-quark-loop plus corrections" at large q'^ , approaching the 
charmonium resonance region. It is the purpose of our paper to address these questions. 

The soft-gluon emission from the charm loop was calculated for the inclusive B — )• 
Xs'j decay in H] and, independently, for the exclusive channel B — )• K*j in Q. The 
main outcome was an effective quark-antiquark-gluon operator whose Wilson coefficient 
is proportional to In i? — t- K*^ the hadronic matrix element of this operator was 

obtained Q using three-point QCD sum rules and, more recently, LCSR Studies of this 
effect in — )• and B — )■ Xsi~^£~^ |8| revealed that, in addition to the lowest-dimension 
local operator, a tower of operators with the derivatives of the gluon field has to be taken 
into account. Explicitly, the soft-gluon momentum was taken into account in for the 
inclusive decay B — Xsi'^£~. 

In this paper we employ operator-product expansion (OPE) near the light-cone for 
the gluon emission from the c-quark loop at <^ 4m^. From this expansion we derive 
an effective nonlocal quark-antiquark-gluon operator. Sandwiching this operator between 
the B and K^*^ states we calculate the hadronic matrix elements employing LCSR with 
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S-meson distribution amplitudes. These sum rules have been used in 11 1 to obtain 
the B — )■ light-meson form factors. The calculation presented below provides an effective 
resummation of the soft-gluon part of the charm-loop effect in i? — t- K^*^i^£~ and (as a 
byproduct) in B ^ K*^. Adding the calculated soft-gluon contribution to the leading- 
order factorizable cc loop , we obtain an estimate of the charm-loop effect in B ^ K^*^£'^i~ , 



valid at < iml. 

The second problem investigated in this paper is the validity of OPE at large q^, 
approaching the charmonium region. This issue is closely related to the status of quark- 
hadron duality for the charm- loop amplitude. Earlier models [12|, adding the contributions 
of ^-resonances on the top of the c-quark loop, involve double counting of quark-gluon 
and hadronic degrees of freedom. More consistent is the use of the hadronic dispersion 
relation written in terms of ip states, as suggested for B — )• Xgl'^i' in [13|. Recently, this 
approach was reconsidered in |14|. In all previous analyses involving dispersion relation, 
nonfactorizable contributions to the charm-loop effect were neglected, and, correspondingly, 
the nonleptonic B — )■ ipK^*^ amplitudes were taken in the factorization approximation. 
In order to adjust these amplitudes to their measured values, additional /c-factors were 
introduced ]l|, [|]. 

In this paper we employ the hadronic dispersion relation in a different way. We fix 
the absolute values of residues of the J/ip- and V'(25')-poles from the experimental data on 
B -)■ J /i;K^*\'4){2S)K^*\ The integral over the spectral density above the open charm 
threshold is parameterized in a form of an effective pole. We then fit the whole dispersion 
relation to the OPE result at <^ 4m^, including the newly calculated nonfactorizable 
contribution. Importantly, we find that this procedure favours a destructive interference 
between the ip{2S) and J/ip terms in the dispersion relation. Our main result for the 
charm-loop effect is obtained in terms of the q^- and process-dependent correction to the 
Wilson coefficient Cg. This correction is directly calculated at small and analytically 
continued to large q^ < fT^'^(2S) OPE-controlled dispersion relation. 

The plan of this paper is as follows. Sect. 2 contains an introductory discussion of 
the OPE for the charm-loop effect. In Sect. 3 we expand the product of the four-quark 
operators with the c-quark e.m. current near the light cone and obtain the non-local 
operator corresponding to the nonfactorizable soft-gluon emission. Sect. 4 contains the 
derivation of LCSR for the B ^ K and B — )• K* matrix elements of this operator. Adding 
them to the known factorizable matrix elements, we obtain the full charm-loop contribution 
to B K'^*k+e-, valid at small q^. In Sect. 5 we present the details of our numerical 
analysis and obtain the results for the charm- loop effect in a form of the correction to Cg . 
Sect. 6 contains our prediction for the charm-loop effect in i? — t- K*'y. In Sect. 7 we perform 
the matching of the OPE result to the hadronic dispersion relation in terms of charmonium 
resonances, predicting the charm- loop effect at q^ up to the open charm threshold. Finally, 
we investigate the infiuence of this effect on the observables of — )• K^*^i'^i^. Sect. 8 
contains the concluding discussion. In Appendix A the effective operators and their Wilson 
coefficients are collected. In Appendix B the definitions of — t- K^*^i'^£'' amplitudes and 
the inputs for form factors are presented. In Appendix C the relevant characteristics of J/ip 
, V(25) states and B K^*\j /^{'4j{2S)) decays are given. In Appendix D the expressions 
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Figure 1: Charm-loop effect in B K^*H'^£ : (a)-the leading- order factorizable contribution; (h) 
nonfactorizale soft-gluon emission, (c),(d)-hard gluon exchange. 

for the coefficients entering the LCSR are presented. 



2. Light-cone dominance of the c-quark loop 



The combined action of the four-quark operators Oi and O2 in (1.1) and the e.m. in- 
teractions of c-quarks and leptons leads to the charm-loop effect depicted in Fig. |l]. The 
contribution of this mechanism to the B — )• K^*^i^i~ decay amplitude can be written as 



-i47raemQc)^VttV,*/-^n^^^'''*'\p,q) , (2.1) 



where Qc = 2/3 is the c-quark electric charge, the lepton current and photon propagator 
are factored out and the hadronic transition matrix element is: 



(2.2) 



Isolating in ( |2.2| ) the T-product of the c-quark e.m. current and the c-quark fields 
entering Oi or O2, one has in both cases a generic expression: 



CiOi{0)+C2O2{0)]}\B{p + q)). 



^xe^''-"r{c(x)7^c(x),CL(0)rcL(0)}, 



(2.3) 
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where denotes a certain combination of Dirac- and colour- matrices. 

At <C 4m^, that is, if c-quarks are highly virtual, the dominant region of integration 
over X in ( |2.3| ) is concentrated near the light-cone = 0. To see that, we follow [11|, where 



the light-cone dominance of the correlation functions containing one heavy-quark field was 
demonstrated. Taking small but time-like, it is convenient to consider the rest frame of 
the virtual photon with q = {qq,^), defining a unit vector w = qj \f(^ = (1,0), w'^ = 1. 
The virtual c-quark fields in ( |2.3| ) are then decomposed into the static part and residual 
4-momentum: pc = rricW + p, so that c(x) = exp(— imcU) • x)hw{x), where the effective 
field hyj{x) contains only the p components. Note that the last redefinition only aims at 
separating the rric scale from the virtual c-quark momentum and making the residual off- 
shell momentum p independent of this scale. Since the virtual c quark 4-momentum is 
Pc = Q — Pc = Q — {mcW + p), the rricW component of the Dirac-conjugated c-field enters 
with the same sign, so that c{x) = exp(— imctf • x)hw{x). After rescaling all c and c fields, 
the operator (|2.3|) transforms to 



C^(9) = I ^i'xe-^(2"^"'-'')-^T{/i^(x)7^/i»(x),/i^i(0)r/i^L(0)} , (2.4) 

where the mc-dependence is now concentrated in the exponent. In the above integral, 
the dominant contribution stems from the region where the exponent does not strongly 
oscillate, that is, where {2mcW — g) • x ~ 1 , yielding 

(x^) ~ l/{2mcw - qf ~ l/(2m, - ^J^f , (2.5) 

(for a more detailed derivation see e.g., |l^). In the region q^ <^ 4m^ the light-cone OPE 
for the product of operators in ( p.4[ ) can be applied. When grows and approaches the 
threshold 4m^, the expansion becomes invalid. Note that due to Lorentz-invariance, the 
estimate (|2.5| ) is valid in any frame, including the rest frame of i?-meson, where both the 
c-quark pair and s quark emitted in the heavy 6-quark decay, are energetic. 

3. Expansion near the light-cone 

The expansion of the operator-product in (2.3) starts with contracting the c-quark fields 



in the free propagators forming a two-point c-quark loop. At leading order, there is no 
difference between the light-cone and local OPE and the hadronic matrix element is: 

4^-^'*')(p,g)]^^^^= (^ + ^2) {K(*\p)\0,{q)\B{p + q)), (3.1) 

where both Oi and O2 contribute and the local operator 

9 

O/^iq) = {q,,qp - 1^9t^p)^^ oi^l, q^)sLYbL ■ (3.2) 



is reduced to the 6 — >• s current. Hence the matrix element ( p.lD is factorized to i? ^ 
K^*^ form factors. Diagrammatically, this contribution is shown in Fig. |l^. The charm- 



loop coefficient function in (3.2) in the adopted operator basis is given by the well-known 
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expression 



-g (2 + !/(<;")) vWPl aictan j , (3.3) 

where = Arn^/q^ > 1, the renormahzation scale is taken as rrih and g{m^,{)) = 

— 8/9 ln(?n,c/mb) — 4/9. Furthermore, a dispersion relation for the coefficient function in 
the variable is valid: 

5 m^,g = g m^,0 + — / — 2^ds. 3.4 

TT J s[s- q^) 

with the spectral density: 



1 . 9 ^ 4 / 4m2 2m^,^, 9, 

-Im,5 me',s =oVl £ 1 + ^ e s-4m2 . 3.5 

vr 9 V s s 



The factorizable amplitude is often called "perturbative" or "short-distance" charm- 
loop effect. There are gluon corrections to this amplitude which do not violate factorization. 
For instance, the perturbative gluon exchanges within the c-quark loop in Fig. [l|a are known 
and can in principle be added in a form of 0{as) corrections to the coefficient function 
g{m'^,q'^). We will neglect them here. The 0{as) corrections to the weak-current vertex are 
implicitly included in the calculation of the hadronic form factors, e.g., in LCSR |17, 

The nonfactorizable contributions to the hadronic matrix element (|2.2D start with 
the one-gluon emission from the charm- loop. Some of the nonfactorizable diagrams with 
hard-gluon exchanges in 0{as) are shown in Fig. |l]c,d. They are included in the NLO 
factorization formula 01, (see also H) for B ^ K<*k+e- 01 B ^ K* 7 in a form of 
the hard-scattering kernels convoluted with the B- and i^*-meson light-cone distribution 
amplitudes (DA's). These perturbative effects can be taken into account separately and 
will not be included in our analysis. Nonfactorizable two-gluon effects in ( |3.2| ) with one 
hard and one soft gluon, or with two hard gluons are also neglected, being presumably very 
small due to extra 0{as) suppression. 

In this paper, we consider the emission of one soft gluon (with low virtuality but 
nonvanishing momentum) from the c-quark loop and evaluate the contribution of this 
mechanism to the B — )■ i^'^*) amplitudes. One of the corresponding diagrams is shown 
in Fig. |l]b, the second one has a gluon emitted from the other c-quark line. 

In terms of the light-cone expansion, the contraction of c-quark fields entering ( |2.2| ) 
yields a gluon-field operator multiplied by a 6 — >• s colour-octet current. Note that only 
the operator Oi contributes at this level. In order to calculate this contribution we use the 
c-quark propagator near the light-cone including the one-gluon term: 

(0 I T{c{xi)c{x2)} \0) = -i [ ^^e-i'=(-i--2).^ + "^'= 



(27r)4 m2 - 



1 
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where G'^^{x) = 5s(A°/2)G""^(a;), u = 1 — u, and the fixed-point gauge for the gluon field 
is used. The above symmetric form can easily be derived from the expression presented, 
e.g., in illl. 

To proceed, we define the light-cone kinematics in the rest-frame of the decaying B- 
meson, introducing the unit vector v = pb/^b = (PkM + q)/^b = (l, 0,0,0), two light- 
cone vectors n±, so that : 

V = 2 (^+ + '^-) ' = = , = 2 , 

Q = {n'.q)'!^ + (n+q)^ , (3.7) 

and choosing = 0. We consider a kinematical situation when the lepton-pair invariant 
mass is small: q'^ = {n+q){n-q) <^ 4m^ < ml , while (vq) = l/2[{n^q) + (n^q)] ~ mi,/2 is 
large, so that one of the components of q in the above expansion dominates, or equivalently, 
q is approximately parallel to one of the light-cone directions. We choose 

2 

q ^+ 

(n^q) mti, (n+g) ~ — ^ (n_g), hence q (n^q) — . (3-8) 

mf, 2 



The propagator (|3.q ) taken between x and involves the gluon field at the point ux, 
(0 < n < 1) which we rewrite via nonlocal differential operator: 

C^^iux) = exp[-iuXr{iV^)] G"^ . (3.9) 

where G"'^ = G"'^(0). Note that due to the fixed-point gauge, the simple derivative can be 
replaced by the covariant derivative acting on the gluon field and ensuring gauge invariance. 

We are interested in the dominant effect of the nonvanishing gluon momenta generated 
by the exponent in ( |3.9D . Decomposing the covariant derivative in the light-cone vectors 

V = (n+P)— + (n_P)^ + V^, (3.10) 

we retain only the n_ component, which corresponds to the gluons emitted antiparallel to 
q, that is, in the same direction as the s-quark in the i?-meson rest frame. We then have 

C'f^iux) ~ exp[-in(n„x)-^^^^^±^]G"^ 

du exp[-m(n_x)w] 6[oj - i!!l±^] G'^^ . (3.11) 



Inserting this expression in the propagator ( |3.6D , and shifting the k integration, we 
obtain for the one- gluon part of the propagator: 



(0 I T{c{x)c{0)} I 0)g = -i I doj I du I -j—^expi-ikx) 







+ mc- uji-u})aai3 + uaap{f. + ruc - uji-uj) _ (m+P) ^^^g 



2(m^ — [k — un^LoY)"^ 
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Using this form of the propagator in ( p. 2]) , we cast the soft-gluon emission part of the 
hadronic matrix element in a form: 



nonfact 



2Ci{K(^*\p)\0^{q)\Bip + q)), 



(3.13) 



where 0^{q) is a convolution of the coefficient function with the nonlocal operator: 



{in+V) 



]Gal3bL 



(3.14) 



In the above, Gajs = -2^apaT 
The coefficient function reads: 



G"'^ , and the derivative acts only on the gluon- field operator. 



I^^pap{q,u}) - ^ 



du 



u%qagpi3 + uqpqagf,i3 - uq Qf^agp/l 



dl{q^,m 



u — u 



-gpaQpli 



I if)}, 



dq 



where 



(3.15) 



dt In 



mi 



t(i-tW 



(3.16) 



is represented in a compact unintegrated form, and we use the notation q = q — uujn^ , 



so that q 



^2 



2uujmb. Here we take into account that w <^ nib, after the hadronic 



matrix element is taken. Note that the neglected components of V in ( 3.1C| ) produce small, 
0{uj/mb) corrections to q'^, hence our approximation is well justified. The spectral density 
of this coefficient function has the following form: 



-Im I i,pc,i3{q,uj) = 

TT 4^^^q^^yq^{q^ - Amj.) Jo 



uqi,qa9pi3 + mpqaOfip 



u + 



[u — u)q 



Ami 



q g^agpp 



-2 Ami) 



(3.17) 



For a gluon field with fixed momentum one restores from ( ^.14 ) and ( 3.15 ) the expres- 
sion obtained in Q from the diagrams with the on-shell gluon emission from the c-quark 
loop. Since only the sum of the two diagrams is ultraviolet-convergent, they have to be 
calculated m. D ^ A. In the resulting expression, contrary to the comment in 1^], we found 
no mass-independent constant term. 

Importantly, the operator ( |3.14| ) has an overall power-suppression factor ~ l/(4m^— g^) 
compensating the gluon field-strength dimension. The B — )• K^*'^ hadronic matrix element 
of this operator is reduced to a specific "nonlocal form factor" 



{K^*\p)\-slY5[u: - ^^^^^]G^MB{P + q)) , 



(3.18) 



which has to be calculated with some nonperturbative method. In fact this matrix element 
resembles a nonforward distribution with different initial and final hadrons. 
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In the local OPE limit, that is, neglecting all derivatives of the gluon field, ( ^.14 ) is 
reduced to: 



with 



and 



= lZpiQ)sLYG^(sbL , (3.19) 



^uLsii) = If^pai3iq, 0) = {qf,qa9pi3 + qpqaQfip - q^ gpagpp)i^^\q^ ,ml) , (3.20) 



^ tii-t) 



At = 0, one easily recognizes in ( 3.ig| ) the quark-gluon operator emerging from the 



charm-loop with soft gluons ^ (see also I, i). The 3-point QCD sum rules in |6| and 
LCSR in 0] were used to evaluate the B — >• i^*7 hadronic matrix element of this operator. 
The same approach was used in ||2^ to estimate the nonfactorizable soft-gluon effect in 
B J/i)K. 

To clarify the difference between the light-cone and local OPE, we return to the non- 



local operator ( 3.14 ) and expand the 5- function in powers of the derivatives acting on the 
gluon field. Integrating the coefficient function over one encounters a tower of local 
operators: 



1 



1^=0 



Sandwiching the n-th operator between B and K^*\ one obtains a contribution of the order 
of (mfeAQ(7£))"/(4m^ — g'^)""*"^. Since itiijAqcd is generally not smaller than m^, there is 
no parametric suppression of the n ^ contributions, as already discussed in More 
recently, the description of the charm-loop contribution in i? — )• Xg'j in terms of nonlocal 
operators was advocated in [p3|| . In what follows, we apply LCSR to evaluate the hadronic 
matrix elements ( p. 18 ), thereby performing an effective resummation of the local operators 



m 



(13:221) 



Expanding the c-quark propagator near the light-cone |20| further, one encounters 
higher-dimensional nonlocal operators with two and more gluon fields, each of them gener- 
ating a tower of operators with calculable, albeit complicated coefficient functions. Simple 
dimensional argument based on the presence of extra gluon fields in i?-meson DA's allows 
one to anticipate an overall suppression of these two and more soft-gluon contributions in 
LCSR by additional powers of the scale l/(4m^ — g^) with respect to the leading one-gluon 
term. 

In exclusive B decays, the soft-gluon effects stemming from the light-cone expansion in 
the framework of LCSR were considered in the analysis of i? — t- vrvr in and B — t- J/^jjK 



in [25 1, where the LCSR with DA's of light mesons were used. In this case as shown in 



1 24], an artificial four- momentum in the vertex of the weak operator has to be introduced. 



In the following section we use a simpler approach where this modification can be avoided. 
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Figure 2: The correlation function used to calculate the B —4- X*-*-* matrix element of the soft-gluon 
emission from the charm loop. The oval denotes the B-meson quark- antiquark-gluon distribution 
amplitudes. The c-quark loop with the emitted gluon generates the nonlocal effective operator Op^. 
The cross indicates the point of gluon emission in the second diagram. 



4. LCSR for hadronic matrix elements 

The B —7- K^*^ transition matrix elements ( |2.2| ) determining the charm-loop effect can now 
be represented as a sum of the hadronic matrix elements of the effective operators and 
O^, containing the factorizable and (soft) nonfactorizable contributions, respectively. 



Let us first consider the B ^ K transition. Adding ( |3.1j ) and (3.13) together, we 
obtain: 



^ + C2 ) {K{p)\0^{q)\B{p + q)) + 2Ci{K{p)\0^{q)\B{p + q)) 



nf^^\p^q) = 

= [{p-q)q,-qh,W^^^\q'), 
where the e.m. current conservation is taken into account and 

n^^^^'Xq') ={^+ C2) A{q') + 2C^A{q') 



(4.1) 



(4.2) 



contains two invariant amplitudes parameterizing the two hadronic matrix elements in 
In the adopted approximation, the amplitude A is factorized: 



A{q' 



^ g{mlq^)f^^{q^), 



327r2 



if one uses ( p.2| ) and the standard definition of the B ^ K form factor: 

{K{p)\-SLlMB{p + q)) = fiK{Q')Pp + ■■■ 



(4.3) 



(4.4) 
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Our main task is to calculate the amplitude A determining the soft-gluon emission 
from the charm-loop. To this end, we employ LCSR with the S-meson DA's and, following 
|1C, 11 1, introduce the correlation function: 

J'S^''\p,q)=i J dSe'P-y{0\T{j!fiy)d^iq)}\B{p + q)), (4.5) 

where = d'^y'^^s is the i^-meson interpolating current and B meson is taken on shell, as 
a HQET state: \B{p + q)) ~ \B{v)). Inserting the full set of states with the kaon quantum 
numbers between the current and in ( |4.5| ) we obtain the hadronic dispersion relation: 



where fx is the kaon decay constant defined as {Q\d'yu^^s\K{p)) = ifxPu and the spec- 
tral density Pufiis,q'^) accumulates excited and continuum states with the kaon quantum 
numbers, located above the threshold Sh- 

Two comments are in order. First, in the approach we are using, hadronic matrix 
elements are related to the correlation function via dispersion relation. Hence, the "full" 
hadronic matrix element B ^ K transition with the soft-gluon emission enters the 
residue of the kaon pole in (^^). In other words, although we have chosen a particular 
correlation function with i?-meson DA's, there is no need to add a contribution where a soft 
gluon emitted from the charm loop enters the final-state kaon "wave function" . Secondly, 
at very large timelike p^ there are also "parasitic" charm-anticharm states contributing 
to the hadronic spectral density in (^), but they are heavily suppressed after the Borel 
transformation in p^. This circumstance allows one to avoid introducing an auxiliary 4- 



momentum in the effective-operator vertex, as suggested in [24|. 



In 1 11 1 the form factor ftjp-iq^) was calculated from the correlation function similar to 



( |4.5| ), where, instead of the complicated effective operator, the 6 — )• s, vector current was in- 
serted. LCSR was obtained including the contributions of two-particle (quark-antiquark) 
and three-particle (quark-antiquark-gluon) i?-meson DA's. Here the leading-order dia- 
grams shown in Fig. || involve only the three-particle DA's of B meson. We calculate these 
diagrams contracting the s-quark fields. The result reduces to the vacuum-to-S matrix 
element of the dGb nonlocal operator. It is decomposed in HQET in four three-particle 
5-meson DA's pi: 



{()\Uy)5[u: - i!!^]G.,(0)6^(0)|^(t;)) (4.7) 
= dXe-'^y^ {1+ /j)S^{v^jr-VrJa)[^A{X,2u;)-^v{X,2u:)] 

Vy J J /3a 

where /b and niB are the S-meson decay constant and mass, respectively. Further details 
can be found in [0] . 
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Equating the correlation function J-yfj,, written in terms of the i?-meson DA's, to its 
hadronic representation (4.6), we perform the Borel transformation in the variable and 



employ quark-hadron duality in the kaon channel to approximate the integral over the 
spectral density py^^, introducing the effective duality threshold Sq. The resulting LCSR 
for the nonfactorizable B ^ K hadronic matrix element of the charm-loop effect reads: 



o-o oo 1 1 2 2 

Aiq^) = ;— — — 7) 7) 7T7 I da Idoj du dt exp ( — ^ , ^„ — — 

87r^fK{ml-ml-q^) J J J J ^ V M2 





1 



C(*^) {q"^, u, a, UJ, t)'^v {mB<y, 2a;) 



ml — t{l — t){q'^ — liJiBUOj) 

+C^^^\q^,u,a,uj,t)XA{mBa, 2uj) - ^[C^^^\q^ ,u,a,uj,t)XA{mBa, 2uj)] 

, n, fj, uj, t)YA{mBcr, 2uj) 

where M"^ is the Borel parameter, do ~ /^"b (^P small corrections), 

XA{X,io)= [^dTXA{T,oo), Ya{X,oo)= [^dTYA{T,Oo). (4.9) 
Jo Jo 

and the coefficients C^'"-* and C^ ' -* are collected in Appendix D. The sum rule is presented 
in ( [1.9D in a simplified form, omitting the s-quark mass and small 0(l/rn-^) corrections. 
In the numerical analysis the complete expression is used. 

Note that the spectral density of the correlation function (4^) in the variable p^, 



in addition to the s-quark pole, contains also "parasitic" contributions which correspond 
to putting simultaneously on-shell the s-quark and the cc-quark lines in the diagram of 
Fig. ^. After Borel transformation these terms are suppressed by ~ exp(— 4m^/M2) and 
we neglect them. As noted above, the corresponding hadronic states are also neglected in 
the dispersion relation. 

In order to assess the accuracy of using the nonlocal effective operator in the 
correlation function we derived the LCSR in an alternative way, inserting the operator 
CiOi + C2O2 and the c-quark e.m. current directly in the correlation function. As a result 
we obtain a more complicated expression for LCSR, which differs from ([4.9| ) only by terms 
suppressed by inverse powers of rric and/or m^, and yields numerically very close results 
for the amplitude A. Finally, the local OPE limit of LCSR (|4.9| ) was also investigated, that 
is, when the soft-gluon momentum is neglected. This limit literally corresponds to putting 
a; — )• in the coefficients C^'"-' and in the denominator of (4.9). The numerical influence of 
this approximation will be discussed in the next section. 

Turning to the calculation of the B — t- K* transition matrix element defined in (|2.2|) 
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we decompose it into the three kinematical structures: 

nlf^^'^Kp^q) = (C2 + ^) {K*{p)\0,iqmp + q))+2C,{K*ip)\d,iq)\B{p + q)) 



3 , 

eM«/37e*"9V'Wi(g') + i[{ml - ml.)e; - {e*-q){2p + q),,]n2{q^) 



^2 



2 —C^P + l)f^ 



— m^* 



n^iq^) , (4.10) 



where e is the polarization vector of the i^*-meson. 

The invariant amphtudes in the above decomposition contain factorizable and nonfac- 
torizable parts stemming from the matrix elements of and O^, respectively: 

■H\^^^'\q^)=(c2 + ^)v^iq^) + 2CiV^iq'') (i = 1,2,3). (4.11) 



The three factorizable amplitudes Vj are easily obtained if one uses (3.2) and the 
standard definition of i? — )• K* vector and axial- vector form factors: 

2{K*ip)\sLjpbL\B{p + q)) = ep^p,e*^qPp^ ^ 



niB + rriK* 
aBK*/ 2\ 

-ie;{mB+mK*)A^''\q^)+i{2p + q)p{e*q) \ ^'^ > +..., (4.12) 

niB + rriK* 

where the form factors multiplying qp do not contribute and are indicated by ellipses. We 
have: 

^2 



Vliq') = 'V ,9imlq')V^^'iq'), 

62TT'^[mB + rriK*) 

V2{q') = -^-^ -n{mlqMf^\q^)^ 

o4:iT^(mB — rriK*) 

,,13, 

b47r^ \mB+mK* niB — rriK* I 

To calculate the nonfactorizable amplitudes Vj(g^) {i = 1,2,3), describing the soft- 
gluon emission we again resort to the LCSR method, introducing the correlation function 

4^,^''*\p,Q)=i I d''ye^Py{0\T{j!f\y)d,{q)}\B{p + q)), (4.14) 

where j^^* = d^y^s is the interpolating current for the K* meson. The hadronic dispersion 
relation for this correlation function reads: 



m'j^t — p y s* s — p'^ 

where fx* is the K* decay constant defined as {Q\j^*\K*{p)) = eymx* fx* and the overline 
denotes the average over K* polarizations. We neglect the total width of K* . This approx- 
imation can in principle be avoided by introducing a Breit-Wigner type parameterization 
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for this resonance. Excited and continuum states with K* quantum numbers above the 
threshold contribute to the spectral density pt^{s, q^). The latter is approximated using 
quark-hadron duality and introducing an effective threshold Sq . 

Furthermore, the K* contribution in ( 4.15| ) is written in terms of invariant amplitudes: 



— 



+qt.qu(v2{q'') - (1 - )V3{q') 



+ ..., (4.16) 



where we only show the kinematical structures that have been used in our analysis. Apply- 
ing the same decomposition to the OPE result for the correlation function T*^ we obtain 
three sum-rule relations which are then used to calculate the hadronic matrix elements 
Vi,2,3('i'^)- The derivation of these LCSR in terms of i?-meson DA's is very similar to the 
one described for B ^ K case. Indeed, the only difference between the two correlation 
functions ( f4.5D and ( [4.14 ) is in the quantum numbers of the s-quark current. The result- 



ing LCSR for Vi{q^) has the same structure as ([4.9|), only the coefficients multiplying the 
S-meson DA's and the overall normalization factors are different. In order not to overload 
this paper, we do not present these expressions here. 

5. Numerical analysis 

We use LCSR with B meson DA's to evaluate not only the nonfactorizable amplitudes A 
and Vi but also the B ^ K and B — )• K* form factors entering the factorizable parts A 
and Vi given by (O) and ( [4.13| ). The sum rules for these form factors and their input are 



taken from ||11[, hence we use the same input for the new LCSR obtained in the previous 



section. In particular, for the three-particle i?-meson DA's the model suggested in [11| is 
taken: 



6^0 



A 



2 



Xa{XM = -^..(2A-a;)e-(^+'^)/'^% 

Ya{\ oj) = -^u;{7u;o - 13A + 3^)6-^^+^^^° . (5.1) 

In this model the parameter uq is equal to the inverse moment A^ of the B meson two- 
particle DA (f>^ and the normalization constant of the three-particle DA's is A|; = 3/2 A^. 



For the inverse moment we use the interval obtained in [27| from QCD sum rule in HQET: 
AB(lGeV) = 460 ± 110 MeV . The scale-dependence of this parameter is neglected. To be 
consistent with the 0{as) accuracy of Xb, the i?-meson decay constant /b = ISOitSO MeV 
obtained from the two-point QCD sum rules in 0{as) is used. Furthermore, the decay 
constants and threshold parameters of K and K* mesons are taken the same as in []ll| : 



Jk = 159.8 ± 1.4 ± 0.44 MeV , /k- = 217 ± 5 MeV , 

= 1.05 GeV^ , s^* = 1.7 GeV^ . (5.2) 
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Figure 3: The nonfactorizable B K amplitude Alji^) of the nonlocal effective operator (solid) 
and its local limit (dash- dotted), plotted as a function of , together with the factorizable amplitude 
A(q'^) (dashed). The soft-gluon contributions are rescaled with a factor of 10. 



For the Borel parameter interval in LCSR we use is = l.OibO.25 GeV^, slightly narrower 
than the interval in [11|. 

In addition, we need the c-quark mass value. Since we are dealing with virtual cc-quark 
loops, it is natural to employ the MS mass, assuming that the normalization scale is about 
2mc- As a default value, we take fhc{2ific) = 1.05 GeV, rescaled from the central value 
of TTic obtained from charmonium sum rule in |28]. To assess the related uncertainty we 
allow the normalization scale to change between 2mc and mc, shifting the c-quark mass 
to mc{fhc) = 1.30 GeV. Finally, m<j(2GeV) = 98 ± 16 MeV is adopted (the average of 
nonlattice determinations |29|). 

In Fig. I our numerical results for the g^-dependence of the dimensionless hadronic 
amplitudes A{q'^) and A{q'^) are displayed for the central values of the input. For conve- 
nience, we normalize both amplitudes to .4(0) = 7.3 x 10~^. The nonfactorizable amplitude 
is about a few percent of the factorizable one and has the opposite sign. Note that both 
A{q'^) and A{q'^) develop imaginary parts above 4m^. In the same figure we plot the nonfac- 
torizable amplitude A{q'^) obtained in the local OPE limit. In this case, the nonfactorizable 
effect changes substantially and diverges approaching to q'^ = 4m^. We come to an impor- 
tant conclusion that the effective resummation the local operators in the framework of the 
light-cone OPE "softens" the gluon correction to the charm loop. 

In the decay amplitude the soft-gluon contribution to the charm-loop effect gets en- 
hanced considerably with respect to the factorizable one by the ratio of the Wilson co- 
efficients 2Ci/(C2 -|- Ci/3) ^> 1. Not surprisingly then, the proportion of the two con- 
tributions depends on the normalization scale of the Wilson coefficients. The scale can 
be fixed more accurately when the perturbative gluon corrections in the hadronic matrix 
elements are taken into account, which is beyond our approximation. For the numerical 
estimates we use the values of the Wilson coefficients given in Appendix A. They are cal- 
culated in leading log approximation at /i = rrih, allowing the scale to vary between O.bmi, 
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Figure 4: The charm-loop effect in B ^ Kt^i"' expressed as a correction to the Wilson coefficient 
Cg (solid), including the nonfactorizable soft-gluon contribution (dashed) with the shaded region 
indicating the estimated uncertainty and the factorizahle contribution (dash- dotted). 
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Figure 5: The same as in Fig. Qfor B^ K*l^l . The results for the correction to the amplitude 
are almost indistinguishable from the ones for AAi shown in the upper panel. 



and 1.5m;,. Note that in our calculation the 6-quark mass value only enters the effective 



Hamiltonian. We adopt the MS value extracted from the bottomonium sum rules in |28| 
TTT-bi^nb) = 4.164 lb 0.025 GeV (conservatively doubling the error). 

It is convenient to express the charm-loop contribution to the B — )• K£^£~ amplitude 
in a form of a (process- and g^-dependent) correction to the Wilson coefficient Cg: 

C79^C9 + ACf'^^^)(g2). (5.3) 
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Substituting and (gj) to and comparing the result with the contribution 



of O9 to the B — )■ Ki'^i amphtude given in Appendix A, we obtain: 



3 f^KiQ') 

(Ci + 3C2) g{mlq') + 2Ci5(^"^'^^^)(g2) , (5.4) 



where the function 



determines the new soft-gluon nonfactorizable part of ACq^'^'^^^^ and represents our main 
result. Here only the ratio of the calculated hadronic matrix elements enter and they both 
are calculated within one and the same LCSR approach. The correction ACg'^'^'^^^^ 
and its factorizable and nonfactorizable parts are plotted in Fig. ^ where the uncertainties 
stemming from our calculation oi'g^'^^'^^^\q^) are indicated. We vary all input parameters 
within their adopted intervals and add individual variations in the quadrature. In Table || 
we display the value and the estimated uncertainties oi'g^'^'^^^^^) at = 1 GeV^, except 
the variation due to uncertainty of which is negligibly small. 

Note that substantial uncertainties are caused by the shift of rric and rather broad 
interval of the inverse moment of the i?-meson DA. 

To parameterize the charm-loop effect for the B — t- K*i'^i^ decay amplitude we use 
its decomposition in the three invariant amplitudes Mi presented in Appendix B and the 
corresponding decompositions ( 4.1[1| ) and ( 4.11 ). As a result the terms proportional to Cg 



in the amplitudes A^i,2,3, have to be modified in the following way: 



function 


g{cc,B^K) 


~{cc,B^K*,Mi) 


g(cc,B^K*,M2) 


~{cc,B^K* Mi) 


cent r. value 


-0.041 


0.26 


0.27 


0.46 


Am, 


+0.014 


-0.08 


-0.09 


-0.15 


Apf2 


+0.00 


-0.04 


-0.04 


-0.07 


-0.001 


+0.07 


+0.08 


+0.12 


Aas 


-0.016 


+0.30 


+0.36 


+0.75 


+0.017 


-0.17 


-0.18 


-0.33 


Atot 


+0.022 


+0.31 


+0.37 


+0.76 


-0.016 


-0.19 


-0.21 


-0.37 



Table 1: The functions determining the soft-gluon correction to Cg, calculated from LCSR, central 
values at = 1 GeV^, and the uncertainties Aa caused by the variations of the input parameters 
(6m^ = +0.25 GeV, Sm- =t°o.ll GeY, 6x, =tll° MeY). 
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where, respectively, 



9 



(cc,B^K\Mi)(^2 



327r^ {niB + mK*)Vi{q^) 
3 qW^^'iq^) 




9 



{cc,B^K',M2)(„2 



647r^ {rriB - rrij^*) V2(g^) 
647r2 [(ms + m|^)V2(g^) 



3 [ 

V3(g2) 



(5.7) 



imB-mK^)Ai^'{q^)\- 



Note that the nonfactorizable contributions to AC( 




are enhanced at small 



q'^ > Arnj with respect to the factorizable ones due to the virtual photon propagator (the 



in i? — )• K*£~^£~ transitions is predicted to be considerably larger than in B ^ Ki~^i~, 
having also a larger uncertainty. Moreover, the nonfactorizable contributions in S — )• K* 
have the same sign as the factorizable ones, which leads to an additional enhancement of 
the charm- loop effect. 

We come to an important conclusion that the calculated charm-loop correction 

remains small, not exceeding (within uncertainties) ~ 5% 
of the Cg value at < < 4.0 GeV^. In B ^ K*£+£- the same ratio may reach (adding 
up the estimated uncertainties) as much as 20% of Cg at 1.0 < < 4.0 GeV^ and is inflated 
at smaller q'^. In both decays, especially in i? — )• K*£'^i~, the nonfactorizable soft-gluon 
part of ACg plays a decisive role. At q'^ > 4.0 GeV^ where OPE for the charm-loop effect 
starts to diverge we will use a phenomenological ansatz for ACg'^^'^^^ "^''\q'^) presented 
below, in Sect. 7. 

6. Charm-loop effect in B ^ K*'j 

We are in a position to predict the charm-loop effect also in i3 — )• K*^ which is simply a 
by-product of our calculation for B — )• K*£~^£~ at q'^ = 0. Importantly, in this case only 
the nonfactorizable soft-gluon emission contributes. With the effective resummation of this 
effect, our result goes beyond the estimates in Q and where the matrix elements of the 
local quark-gluon operator were calculated. 

To specify the normalization, in Appendix B we present the dominant contribution to 
the B — )■ K*j amplitude, due to the operator O7. The charm- loop effect is included in this 
amplitude by the following additions to the Wilson coefficient C^^-^ multiplying the first 
and second kinemetical structures in this amplitude, respectively: 





(6.1) 
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where 

..^{ccB^K'^). _ 327r2 CiVi(O) 

^ ' 3 (mb + m,)Tf^*(0) ' 

64^2 CiV2(0) ,^ 

The numerical analysis reveals that the corrections to both amplitudes are approxi- 
mately equal: 

^^^(c-c,B^X*7)] ^ ^ ^^^{ccB^K*^)^ ^ ^ (-1.2l?|) X 10-2 , (6.3) 

amounting up to 8 % of the coefficient Cj^^ . Our estimate yields larger magnitude (also 
with larger uncertainty) than predicted in Q where the local OPE and LCSR with K*- 
meson DA's were used. Normalizing the hadronic matrix elements calculated in Q as in 
( ]6.2| ) and using their input, we obtain 

[Acf '^^^*^)]f ^ = (-0.39 ± 0.3) X 10-2 , 

j^^{c-c,B-^/<*7)]p ^ (_o g5 ^ 0.57) X 10-2 . (6.4) 



Note that if we also use the local OPE limit, both amplitudes in (|6.3D have about 40% 
larger magnitudes, approaching the estimates obtained from the three-point QCD sum 
rules in 

7. Accessing large with dispersion relation 

Returning to the hadronic matrix element (^^) we use its analyticity in the variable q^, 
expressed in the form of dispersion relation. Considering first the B ^ K case, we obtain 



the corresponding spectral density, inserting in (2.2) the full set of hadronic states with 
the quantum numbers of J/ip between the c-quark current and the four-quark operator. 
Introducing the hadronic matrix elements 

mYc\^{q)) =e';m^U, (7.1) 

{K{p)i,{q)\CiOi + C202\B{p + q)) = (e^ • p)m^AB^K , (7.2) 

where ip = J/ip,ip{2S) and is the polarization vector, we obtain the dispersion relation 
for the invariant amplitude: 



oo 



+ / ds- '^'^ 



Ami - 9^ - ^e) 



(7.3) 



where p{s) is the spectral density of ^-resonances and continuum cc-states located above 
the open charm threshold, s > Arn?j^. In ( [7.3| ) we include the small total widths of J/ip 
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and 'tp{2S), but neglect the complex phase everywhere beyond the immediate vicinity of 
the resonances. One subtraction at = takes into account that the spectral density 
( ]3.5| ) of the factorizable loop contained in T-L^^^^\q^) tends to a constant at s— s-oo. Note 
that the soft-gluon contribution has a convergent dispersion integral as follows from the 
spectral density ( p. 17 ). The parameters of V'-poles in (O) are presented in Appendix C, 



including the decay constants and the absolute values of the invariant amplitudes AstpK 
calculated from the measured B — )• ipK widths. 

It is well known that the naive factorization approximation for these amplitudes: 
Abi/jK = {C2 + Ci/3)f^f^j^{m'^) is not consistent with the experimental data, indicating 
sizeable nonfactorizable corrections. Moreover, both amplitudes A^j/^k ^iid Ab^(^2S)k ^re 
expected to have complex phases generated by the strong final-state interactions in these 
nonleptonic decays ^. Note that the diagrams with perturbative gluons (e.g., the 0{as) 
diagram in Fig. lc,d) generate complex phases in 'H^^'^{q'^). One can speculate that they 
are dual to the final-state interaction phases in A^ijjKi in a certain analogy with the QCD 
factorization approach |3^] . The diagrams with perturbative gluons are not included in our 
calculation of 'H^^^{q'^), hence we also neglect complex phases in the hadronic part of the 
dispersion relation. On the other hand, we notice that after taking the soft-gluon emission 
into account, the amplitude 'H^^^^\q^), is not positive-definite anymore. Indeed, the 
nonfactorizable part A of this amplitude obtained from LCSR has a negative sign with 
respect to the factorizable part A. Therefore, we relax the positivity condition for the 
spectral density in ( [7. 3D , allowing, e.g., different signs between the residues of ^-poles ^. 

in. 



This makes (|7.3D essentially different from the dispersion relations used earlier in [13, 

The most complicated part of the dispersion relation ( [7.3[ ) is the spectral density p{s) 
above the open-charm threshold. It includes an interplay of broad charmonium resonances 
and continuum states. Since we are working beyond the factorization approximation, we do 
not attempt to describe this spectral function as a sum over higher resonances and a tail 



determined by quark-hadron duality with the c-quark loop as in [13] . In future p{s) can be 
partially determined and/or constrained up to s = {ms — m^)^, employing the experimen- 
tal data on the nonleptonic decays B — )- ipK, {ip = ^(3770), ...) and B — > DDK, D*DK, .... 
Note also, that at q^ approaching the upper threshold {niB —rriK)'^, the dispersion relation 
will also be infiuenced by the singularities in p{s) related to the intermediate hs states. 
Hence we refrain from using the quark-hadron duality approximation for the integral over 



p{s) in (7.3). The simplest possible choice is to model this integral by an effective pole: 



L 



00 ^(g) a^B^^) 



,2 s[s - q^) mi - q^ 



The next step is to match the dispersion relation to the function 'H^^^^\q'^) obtained 
from OPE and LCSR at q'^ <C 4m^. More specifically, we use as an input the l.h.s. of 



^ These phases originate from the discontinuities of the B — > ipK amplitudes in the variable (p + q)^ and 
are not related to the analytical properties in the variable g^. 



^The sign of the lowest J/ip contribution with respect to the calculated amplitude in l.h.s. of (7.3) 
is simply fixed from the higher derivatives of the dispersion relation over at small q^ where all other 
contributions in r.h.s. are power suppressed. 
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dfj) calculated at -Ami < Q'^ < 2 GeV^, so that the subtraction constant H^^^I^O) is 



also determined. Provided the absolute values of t/'-poles in (7.3) are fixed by experimental 
data this matching allows one to fit the parameters of the effective pole. Importantly, 
the fit clearly favours a negative relative sign of J/ip and tp{2S) contributions, and yields 
for the central values of the input 



4.06 GeV, a(^^^) = 0.06 x 10"^ . (7.5) 



The fitted residue of the effective pole turns out to be much smaller than residues of J/ip 
and ip{2S) in (p|): 

''"^^^ ^^^^^ 



To investigate the sensitivity of the dispersion relation to the adopted ansatz ([7.4|), we also 
used, as an alternative option, the conformal mapping — t- z and the z-expansion for 
the integral ( |7.4| ), making use of the fact that it has no singularities at < Arnj^. The 
resulting dispersion representation for T-L^^^^\q'^) obtained after fitting the coefficients 
of z expansion numerically differs very little from the one with the effective pole, again 
favoring the sign pattern as in (^^), hence we adopt the latter as a default model. 

Finally, we express the dispersion representation of the charm-loop effect in terms of 
the correction ACg^^'^^^\q'^) defined as in (|5.4D . Our numerical prediction is plotted in 
Fig. ^ up to q^ = rn^^2s) ■ At < 4 GeV^ it coincides with the calculated result shown in 
Fig. ^. The dashed region in Fig. ^ indicates now all uncertainties, in particular the one 
which is not related with our calculation and corresponds to the substantial scale- variation 
of the combination of Wilson coefficients in the factorizable part. Note that at > 4 
GeV^, where we rely on the dispersion relation, it is not possible to split ACg^'^'^^^^ into 
factorizable and nonfactorizable parts. The predicted charm- loop effect in i? — )• Ki^i~ 
is numerically unimportant at least up to ~ 5 — 6 GeV^. Within uncertainties it is 
even consistent with zero, due to possible cancellation of factorizable and nonfactorizable 
contributions to ACg{q'^)^'^'^'^^^^ at certain combinations of the input parameters. 

Numerically, we find that the correction to Cg is well reproduced by the following 
simple parameterization valid at < g^ < 9 GeV^ 

AC7f'^^^)(0)+ri^ 



Acf '^-^)(,^) = ^^^j^, ^ , (7.7) 



with the calculated 



and fitted 



ACf'''-''')(0) = 0.17l[!;;^ (7.8) 



r, = OMt'oll ^2 = 1.021^1- (7.9) 



This representation can be used in the phenomenological analysis of i? — t- K£'^£~. 

The dispersion approach described above is extended to the charm-loop effect in 
B — )• K*i^£~ , in which case the hadronic matrix element ( p. 2D consists of three invariant 
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Figure 6: T/ie charm loop contribution to the Wilson coefficient Cg for Bq Kl^l^ at below the 
open charm threshold, obtained from the dispersion relation fitted to the OPE result at ^ 4m^. 
The central values are denoted by dashed line, shaded area indicates the estimated uncertainties. 
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Figure 7: The same as in Fig. I for Bo KH+1- 



amplitudes \q'^)j ^ = 1; 2, 3. For the first two ones we employ dispersion relations 
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with two subtractions (taking into account one extra power of in the factorizable part): 



(B^K") 



He? 



dq- 



(0) 



(7.10) 



E 



(z = l,2) 



whereas for the third amphtude a combination of two dispersion relations has to be used 
yielding: 



n 



fi>^BlpK* 



V.=J/V>,V'(2S) 



+ hsiq^ 



(7.11) 



In the above, A^^^x* invariant amplitudes determining B — t- ipK* nonleptonic de- 

cays. They can be expressed via transversity amplitudes (see Appendix C for the definitions 
of the latter): 



A 



(1) 

BipK* 



A 



(3) 

B'il)K* 



I'lB^BipK* 

2mK*m^A%^j^, 



.(2) 

^BijjK* 



-A 



BipK* 



V2{ 



m 



B 



m 



K*J 



A 



m 



B 



m 



K* 



A 



171% X 



+ 



B'^BipK* 



.(7.12) 



The amplitudes A^^^J and their relative phases are extracted using the kinematical anal- 
ysis of i? — ;> Jj'^K* and B ^(2S)K* decays (see, e.g.,|Q), together with the latest data 
on the angular distributions The integrals over the spectral density of higher states 
denoted as hi{q^) in ( [7.11 ) and (7.11) are parameterized with the help of effective poles: 



hi{q^ 



a- 



(B^K") 



(7.13) 



After that the dispersion relations are fitted to the calculated 'H\^~^^ at q^ <^ 4m^. 
Without going into further details, we only mention that in this case the pattern of relative 
signs and hierarchy of contributions is very similar to the B ^ K case. In particular, the 
position of the effective pole coincides within small uncertainties with the one in (|7.5| ) . 

With the help of the dispersion relation we finally obtain the corrections ACg ' ' ' 
presented in Fig. |^ 

For the phenomenologically interesting region 1.0 < q^ < 9.0 GeV^ we suggest the 
following numerical parameterizations of these corrections: 



AC, 



{cc,B^K*Mi) ^^2 



^{B^K*,Mi) ( ^ 



1 + r. 



{B^K*,Mi)q2-g-2 



(7.14) 



where the calculated values of AC, 



{cc,B^K*,Mi) 



"i^2^^ '^'^ collected in Table 



9 



if 



IGeV^) and the fitted values of 
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Figure 8: The differential width of Bq K^^^^ , including the charm-loop effect calculated with 
the central values of input (solid, the shaded area indicates estimated uncertainties) and without 
this effect (dashed). 



parameters 


value at g2 = 1.0 GeV^ 


'1 


{B^K*,Mr) 

'2 


^^(cc,B^K*,Mi) 
^^{cc,B^K* ,M2) 
^^icc,B^K* ,M3) 


' ^-0.37 
''-'-0.41 
J^--L^-0.70 


"-'•J^'-'-o.oo 

'-'•^^-0.00 

'-'•'-'O-o.io 


1 iq + 0.00 

^•-•^"J-o.oi 

1 19+0.00 

^•-•^^-o.oi 
1 n'^+o"^ 

^^•'-'^-0.04 



Table 2: Parameters of the ansatz (7.14) for ACg 



¥K*,Mi) 



(q2) m^id at 1.0 <q^ < 9.0 GeV^ 



In the region between J/'4> and il>{2S) our results for ACg'^'^"^^* '(^^) provide at most 
crude estimates. Nevertheless, this region is interesting because the predicted destructive 
interference between the J/ip and ip{2S) poles manifests itself in the form of a charac- 
teristic maximum located in the middle. In case of the constructive interference (which 
is demanding an unnaturally big, destructively interfering contribution of higher ^/^-states 
in order to satisfy the dispersion relation) the maximum is replaced by a monotonously 
increasing curve. 

At q'^ > "1-^(25)' dispersion relation becomes too complicated to be treated by any 
simple model, and the estimate of the charm-loop effect remains an open problem. From 
what we discussed above, it is obvious, that neither the approximation { cc-loop © gluon 
corrections}, nor a simple sum over ip resonances can provide an adequate description of 
this effect in S ^ K'^*k+£-. 

We are now in a position to investigate the impact of the predicted charm-loop effect 
on the observables in B ^ K^*^i'^£~ . With the decay amplitudes defined in Appendix B 
the differential widths are calculated adding the charm- loop corrections to Cg. Since we 
are only interested in this effect, the other small nonfactorizable contributions (e.g.. the 
loops due to the quark-penguin operators or u-quark loops) are not taken into account. 
Remember that our analysis also does not include perturbative nonfactorizable effects, 
hence our calculated widths are strictly speaking not yet the complete predictions to be 
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Figure 9: The differential width of So 
are the same as in Fig. 8. 



K* ji^ ^ normalized at ~ 1.0 GeV^, The notations 



compared with the data. Furthermore, for the calculation of the B — t- Ki^i" differential 
width we use the B ^ K vector and tensor form factors obtained from more accurate LCSR 
with kaon DA's |jl^, |3^ (see Appendix B). The result is plotted in Fig. ^ As expected, 
the charm- loop effect becomes essential only if one approaches the J/ij: resonance region. 
The same effect is significantly more pronounced in the differential width of S — >• 
where in each Mj-part of the decay amplitude one has to replace Cg by Cg + 

{cc,B^K* ,Mi) 



^^•^L.,a—ri\ ^ The result is presented in Fig. y where we normalize the width to its 
value at = 1 GeV^ in order to diminish the large uncertainties contributing to the width 
by — 7- K* form factors. The latter are calculated from LCSR with B meson DA's (see 
Appendix B). 

Finally, we consider the forward-backward asymmetry in B ^ K*£~^£~ (for definition 
see e.g. |15] ). In our notations, the zero-point of this asymmetry is determined by the 



following equation 



r^eff ^Bimb + rus 

^7 2 



1 + 



+ 1 



rriB 



rpBK' 

^2 



j^BK' 



1 



{cc,B^K*,Mi), 



9 



(%') + AC. 



{cc,B^K'' ,M2) ( 2 
9 



{it) 



0. 



(7.15) 



The numerical result for the asymmetry is plotted in Fig. 10 
equation and taking into account the uncertainties we obtain 



Solving the above 



ll 



2.9tO-2GeV^ 



'-0.3 



(7.16) 
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Figure 10: The forward-backward asymmetry for Bq — > K* decay. The notations are the 
same as in Figs 7,8. 

This is close to the prediction of [|l5| and the one (at the leading-order) of Q. Since 
the above interval at least marginally belongs to the validity region of QCD calculation, 
the estimate ( [7.16| ) is practically independent of the dispersion-relation ansatz. Hence, we 
can clarify the influence of the nonfactorizable contribution on the position of the zero- 
point. Switching this contribution off, we observe a small shift of the central value towards 
(7q = 3.2 GeV^. Note that, according to Q, the addition of perturbative NLO corrections, 
causes a significant shift of upwards. 

8. Conclusion and outlook 

In this paper we went beyond the perturbative approximation for the charm-loop effect 
in — ?• K^*^l'^l^ . In addition to the leading-order c-quark loop, the long-distance soft- 
gluon emission which violates factorization were taken into account. Employing the OPE 
near the light-cone, we derived a nonlocal effective operator responsible for the soft-gluon 
emission from the intermediate c-quarks. The B — )• /C*^*) hadronic matrix element of this 
operator has a structure of nonforward (flavour changing) parton distribution, and has 
to be obtained from nonperturbative QCD. We employed LCSR with B- meson DA's to 
calculate these matrix elements. The same method and input was used to obtain the 
B —7- K'^*'^ form factors entering the factorizable contribution. In the decay amplitudes, the 
new soft-gluon contribution becomes numerically important, being enhanced by its Wilson 
coefficient and by ~ ^/(f' for B — )• K*i^i~ with respect to the factorizable part. The 
characteristic power suppression of the soft-gluon contribution is ~ l/(4m^ — q"^), signaling 
that the approximation "perturbative loop © soft-gluon corrections" is only applicable at 
<C 4m^. At q^ approaching the cc-threshold, multiple soft-gluon emission operators 
have to be included, new hadronic matrix elements proliferate and one eventually looses 
control over OPE. A clear footprints of these effects are the observed large violations of 
factorization in i? — t- charmonium -|- K^*^ decays. 
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LCSR with S-meson DA's used in this paper are not yet sufficiently accurate. There 
are unaccounted ~ l/nib corrections to the HQET correlation function, the gluon radiative 
corrections are not taken into account and the parameters of i3-meson DA's, e.g., the 
inverse moment, still have large uncertainties. The imperfection of sum rules eventually 
converts into relatively large theoretical errors of our calculation. Nevertheless, since we 
are investigating a small effect, at least at low q'^, the achieved accuracy is reasonable. 
We believe that in future also other methods, first of all, lattice QCD can be used to 
calculate the B — )■ K^*^ hadronic matrix elements emerging from the light-cone OPE. Also 
the conventional LCSR with light-mesons DA's and B -meson interpolating current can be 
employed for this purpose, with necessary modifications [^] (an artificial 4-momentum in 
the vertex of the effective operator). 

In this paper we also applied the OPE-constrained dispersion relations for the B — )• 
K^*^ hadronic amplitudes with charm loop in order to to access the large region, including 



charmonium resonances. Our approach goes beyond the earlier ansatz [13, ITH] in terms of 



dispersion relation, since we include nonfactorizable corrections and use the QCD result 
at small and spacelike q"^ to constrain the integral over higher states. This analysis clearly 
indicates a nontrivial interference between the contributions of J/?/; and ip{2S) states. 

One conclusion of our study sounds rather pessimistic. In our opinion, it is difficult, if 
not impossible at all to make a reliable prediction for the charm-loop effect above 'tp{2S), 
based on QCD. Although the actual effect could be small in this region, it will depend on 
the interference of many charmonium states and cannot be reliably constrained by OPE. 
In principle, new accurate data on the branching fractions B — )• ^pK^*^ for states heavier 
than ip{2S) and eventually also on i? — DDK^*^ can be used to saturate the spectral 
density of the charm loop amplitude at large g^. 

In this paper we considered only one, albeit important effect of four-quark operators 
with c quarks in B — t- K^*^i~^i~. Similar effects of quark loops with different fiavours 
stemming from CKM suppressed and/or quark-penguin operators can be taken into account 
along the same lines. Note that in the case of light-quark loops, e.g., the u-quark loops 
important in i? — ?■ p{7r)i~^i~ , one has to expand at large spacelike q^. A separate interesting 
question deserving future study is the potential role of soft gluons in the weak annihilation 
contributions in rare semileptonic B decays. 
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Appendix A: Effective Hamiltonian 

For convenience, here we list all those operators Oj and their Wilson coefficients Cj, entering 
H^f f in (|l.lD which are used throughout this paper omitting Os-e with very small Wilson 
coefficients and Osg ■ 

where the notations gi(/j) = ^—^"^^q and L{R) = ^~Mi5 used. We use the standard 
conventions for the operators Oj, except the labelling of Oi and O2 is interchanged. 

The sign convention for Oj-y and Ogg corresponds to the covariant derivative i-D^ = 
id^ + eQfAfj^ + gT^'A^j^, where Q/ is the fermion charge. In addition, the convention for the 
Levi-Civita tensor adopted in this work is Tr{7'^7''7^7'^7^} = 4,ie^^'^f^ , e°-^^^ = —1. 

We use the Wilson coefficients Ci calculated in the leading approximation. They are 
given in Table ^. 



Appendix B: Decay amplitudes and 
form factors 

1. B ^ Kl+l- 

The dominant contributions to the ampli- 
tude B —7- Ki'^l^ stem from the hadronic 
matrix elements of the operators Oj-y, O9 and 
Oio: 



fi (GeV) 


O.bmb 


mb 


l.bmb 


Ci 


1.180 


1.117 


1.090 


C2 


-0.380 


-0.267 


-0.214 




-0.363 


-0.319 


-0.298 


C9 


4.435 


4.228 


4.034 


Clo 


-4.410 


-4.410 


-4.410 



Table 3: Numerical values of the Wilson coef- 
ficients at three different scales. 



A{B K£+t 



Gp Oien 



^/2 vr 
2{mb + ms) ^eff rT 



+ 



m-B + rriK 



i.l) 



where we omit the contributions yielding the lepton mass (the latter are included in the 
numerical analysis). The tensor B ^ K form factor is defined as 



{K{p)\-sa^pqPb\B{p + q)) 



q^{2p^ + 9p) - (m| - m\)q^ 



ruB + niK 



i.2) 
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2. B^K*l+i- 

For B — >■ K*£'^£~ decay amplitude we use the following expression: 



A{B K*rr 



-ie*^M2{q'^) + lie* ■ q)Pt,M^{q^) 

+h^l^l [w<^e*'^9Wi(9') - ielM2{q^) + i{e* ■ q)p^^^{q^) 



where 



Mi{q^) = 



ffmb + ms BK' (2 



ins + rriK* 



' rpJ: 

.2 "^1 



{1^ 



Msiq^ 



2Cc 



'rriB + rriK* 



and 



Af^'iq') 



2Cio 



10- 

The tensor B K* form factors are defined as 



m o — m 



-(2p + g)^ 



withTf^*(o) = r2^^*(o). 



3. B ^ K*j 

The amplitude of B — > K*7 decay is 

Gf e 



A{B ^ i^*7) = -4^^tbVZC'/^{ - {m, + ms)e 



V2 27r2 

+z(m, - m,) [(e* • e;){p ■ q) - (e* • g)(e; • p)] }rf (0) , 



where e* is the polarization vector of the photon. 



4. B ^ form factors 

The B — )■ K^*^ form factors used in the factorizable parts of the charm-loop effect cal- 
culation are substituted by the corresponding LCSR with B meson DA's. For the dif- 
ferential widths we need an explicit numerical parameterization of these and also tensor 
form factors. We parameterize the g^-dependence of all B — t- K^*^ form factors with the 
z-parameterization similar to the one suggested in [35|, with one slope parameter: 

^(9') = 1 ^2 j 1 + ^1 ,to)-z{0, to) + ^ [z{q^tof -z{0,tof])\, (8.8) 

i q /"i^^(jp) L V ^ /J 

and 



\lT^-q^ - y/T+ - To 



To = T+ - - T-y/r^ . (8.9) 

The pole corresponds to the sh- resonance with appropriate . The B ^ K form factors 
are calculated from more accurate LCSR with kaon DA's |4| at < 12 GeV^ including 
the spacelike region, and the result is fitted to the above parameterization. The B — t- K* 
form factors are calculated from LCSR with B meson DA's in the same region. For 
the two tensor form factors T2^3 the results are new. The form factor values at = and 
the fitted slope parameters are presented in Table ^ where the masses of the Bg states 
are g]: mB,(o-) = 5.366 GeV, mB-^^i-) = 5.412 GeV and mB*(i+) = 5.829 GeV. Within 
uncertainties, the form factors obtaned from the two different LCSR agree with each other 
(seeQ). 

Appendix C: Parameters of ^ = {J /if) ^ if) {2 S)} and B — > ^K^*^ amplitudes 

The necessary data on the two lowest charmonium levels are collected in Table ^ where 
also the absolute values of the B — )• ■i/'X decay amplitudes are given, calculated from: 

^A,,.,^J^lEM^m\"\ (8,0, 
\T(B)Gj.|l/„|2|V;t|2m|A||,J 

where XstpK = A(l, mj^ /m'^,m'^/m'^), and A(o, b, c) = a?' + iP' + (? — 2ab — 2ac — 2bc. 
For B — )• tl^K* we use the decomposition in transversity amplitudes: 

{K*ip)^iq)\ (CiOi + C2O2) \Bip + q)) 

^/2e„e^i^6°/^^>,g.4%^ + 9'-'^^ (8.11) 
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where Xb^k* = X{\,m\t /m\,m'^/m\), so that the decay width is 



i=0,||,± 



The polarization fractions fi{i = 0, ||, _L) are defined as follows: 

/. = — 1 ■ (8-13) 

The transversity amplitudes in this decay can be determined from the data on the 
branching fractions and polarization fractions collected in Table H- In addition, one can 
determine the relative sign between Aq and Ah. 



Appendix D: The coefficients in the LCSR 

The resulting expressions for the coefficients C^^^q"^ ,u,a,u},t) entering the sum rule ([4.S|), 
F = ^v, ^AV, ■■■ are presented here in a form of the quadratic polynomials in the variable 
a;, 




form factor 


pi 
BK 


w(0) 


b\ 




input 












at g2 < 12 GeV^ 




0.34 


+0.05 
-0.02 


9 1+0.9 
~^-J^-1.6 






fO 

JBK 


0.34 


+0.05 
-0.02 


*-3-0.9 


no pole 


LCSR 


JBK 


0.39 


+0.05 
-0.03 


9 9+1.0 
^•^-2.00 




with K OA's 


yBK' 


0.36 


+0.23 
-0.12 


-4.«_o.4 


B*s{l-) 




Af^' 


0.25 


+0.16 
-0.10 


n 34+0.86 
'-'•"J^-0.80 


Bs{l^) 






0.23 


+0.19 
-0.10 


-0.85t?:| 


Bs{l+) 


LCSR 




0.29 


+0.10 
-0.07 


-18-2iko 


Bs{0-) 


with B DA's 


r^BK* 


0.31 


+0.18 
-0.10 


_4 6+0.81 

^•"-•-0.41 


B*s{l-) 




rnBK* 


0.31 


+0.18 
-0.10 


q 9 + 2.1 
<J-^_2.2 


Bs{l+) 




TBK* 


0.22 


+0.17 
-0.10 


-10.3+3:^ 


Bs{l^) 





Table 4: The B — >■ K^*^ form factors from LCSR and their z-parameterization. 
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The nonvanishing coefficients of these polynomials are: 



C^^^> = (1 - t)t [2q\{2a - 5)n + 2) - 2m|(a - l)\u - 1)] 

-^q\2a - 3){2u - I) , 

Cj*^) = 4(1 - t)tu [m|(a(2(fT - A)u - a + 5) + 7u - 4) + 2q^{u - 1)] 

+m|(2cj - 3)u{2u - 1) , 
C^*^) = 8m|(l - t)tu^{a{Au - 3) - 3u + 3) ; (8.14) 



= (1 - t)t[ml{a - l)2(n - 2) - g2(6cTU + u - 2)] + q^a{2u - 1) , 
fji^Av) ^ _ 2{l -t)tu f ^2^^^ _ r (2^^^ _ 1) _ 7^ + 4) + ^ _ 2] 

C7 — 1 \ 

+g2(2a + u-2) ) +2m|o-(l-2u)u; (8.15) 



[XA) 



(Xa) 



(1 - t)t [m|(o- - lf{u - 2) + g2((5 _ 6cj)ii - 2)] 
mB{a - 1) 
q^{2a - l){2u - 1) 
2mB{fT — 1) ' 
4(1 - t)tu [m% {{a - 1)2 - {{a - 4)o- + 2)u) + ^^(i _ u)] 



+ 



niBicr - 1) 



+ 



mBu{2u - 1)(1 - 2a) 
a-1 



(8.16) 



state 






Mass (MeV) 


3096.916 ±0.011 


3686.09 ± 0.04 


Ttot (keV) 


93.2 ±2.1 


317 ±9 


Til (keV) 


5.55 ±0.14 ±0.02 


2.38 ±0.04 


U (MeV) 


416+^ 


297+^ 


5i?(S° ^ J/V^/fO) 


8.63 ± 0.35 X 10"^ 


6.55 ±0.66 X 10""^ 


I^B^^I(MeV) 


3i;? 


411^ 


BR{B° J/iIjK*^) 


13.3 ±0.7 X 10"^ 


7.10 ±0.62 X 10^"^ 


fo 


0.556 ±0.009 ±0.010 


0.48 ±0.05 ±0.02 


f\\ 


0.211 ±0.010 ±0.006 


0.22 ±0.06 ±0.02 


fi. 


0.233 ± 0.010 ± 0.005 


0.30 ± 0.06 ± 0.02 


l^oKGeV^) 


n 97+0.00 
U-^'-O.Ol 


'-'•^^-0.02 


|A|||(GeV3) 
|A±|(GeV3) 


U-J^D-O.OO 


1 5+°-°2 




17+0-02 

'-'•^-'-0.02 



Table 5: The characteristics of J/^ and ip{2S) taken from |3^ . 
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^ 2q\l - t)tu[ml{a - if + q\2a - 1)] 
m%{a-l) 
q^{2u - 1) [m|(a - 1)^ + q^{2a - 1)] 
2m|(cr - 1) 
^ n(2n-l) [m|(c7-l)^ + g^(2c7-l)] 
^ mB(cr-l) 

4(1 - [^|(^ _ 1)2 + q2^2a - 1)] 



.17) 



mB{o - 1) 

C^^-^) = -2mB((j - 1)(1 - t)t{u - 2) , 
= 8mB(c7 + 1)(1 - t)t{u - l)u , 

C7f-) = 8ms(l - t)tu' {j^^^ + 4) ; (8.18) 



C^^-^) = —(a- 1)(1 - i)i(m|((7 - l)\u - 2) 

TUB 

2<n nx\ 2g2(a- 1)(7(2U- 1) 
-g^(6aM + M - 2)) + '—^ '- 



iriB 



C[^' = (1 - t)tu[ml{a - l){a{a{u - I) - Au + 3) + u - 2) 

rriB 

+q^{au + a + u-2)] - iniBicr - l)a , 

C^^^) = 1 -(l-t)tu^\ml(a-lf(2a(u-l)-u + 2) 

mB{a-l) 

+g2(2(T + u-2)] . (8.19) 
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